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O : ABSTRACT 

o 

■ We have observed the quadruply lensed quasar IRXS Jl 131-1231 with the 

pLn ' integral field spectrograph mode of the Kyoto Tridimensional Spectrograph II 

mounted on the Subaru telescope. Its field of view has covered simultaneously the 
three brighter lensed images A, B, and C, which are known to exhibit anomalous 
flux ratios in their continuum emission. We have found that the [OIII] line flux 
^ ' ratios among these lensed images are consistent with those predicted by smooth- 

0> ■ lens models. The absence of both microlensing and millilensing effects on this 

^ ■ [OIII] narrow line region sets important limits on the mass of any substructures 

^ ' along the line of sight, which is expressed as Me < lO^M© for the mass inside 

O ■ an Einstein radius. In contrast, the H/3 line emission, which originates from the 

• broad line region, shows an anomaly in the flux ratio between images B and C, 
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i.e., a factor two smaller C/B ratio than predicted by smooth-lens models. The 
ratio of A/B in the H/3 line is well reproduced. We show that the anomalous C/B 
ratio for the H/3 line is caused most likely by micro/milli-lensing of image C. This 
is because other effects, such as the differential dust extinction and/or arrival time 
difference between images B and C, or the simultaneous lensing of another pair 
of images A and B, are all unlikely. In addition, we have found that the broad 
H/3 line of image A shows a slight asymmetry in its profile compared with those 
in the other images, which suggests the presence of a small microlensing effect 
on this line emitting region of image A. 

Subject headings: dark matter — gravitational lensing — quasars: emission lines 
— quasars: individual (IRXS J1131-1231) 



INTRODUCTION 



Gravitational lensing is an important probe for the r nass distribution of len sing galaxies 
as well as for the internal structure of lensed sources (e.g.. lKochanek et al.ll2006l ). Quadruply 



lensed quasars are particularly useful in this regard since the positions and flux ratios of the 
lensed images provide strong constraints on lens models, and provide far more observational 
information than doubly imaged systems. For instance, recently revealed quads with "anoma- 



Metcal: 



Keeton 



ures in the lensing galaxy (e.g.. 


Mao & Schneider 


1998; 


Dalai & Kochanek 


2002; 


Schechter & Wambsganss 


I2OO2 



2OO3I ). Whether or not these substructures correspond to cold dark matter (CDM) 



subhalos (milli- lensing) , stellar populations (micro- lensi ng), or differential reddening by dust 



remains an issue, and depends on the lens system (e.g.. lChiba et al.ll2005l ). It is of particular 
interest to assess whether there is a high level of substructures or subhalo s in galaxy-sized 
halos, as predicted in CDM models (IKlypin et al.lll999l : iMoore et al.lll999l ). 



A technique to distinguish the origin of anomalous flux ratios was proposed by lMoustakas &: Metcall 



(120031 ). utilizing spectroscopic information of lensed quasars. If we take into account a finite 
(i.e., not point-like) source size of a quasar heart, which consists of a small continuum emit- 
ting region (CR) with typical size i?s of ~ 10~^ pc, broad line region (BLR) with i?s ~ 1 pc, 
and narrow line region (NLR) with i?s ^ 100 pc, then the lensing magnification will be 
functions of these source sizes and the mass of any substructure inside its Einstein angle. 
Me (i.e., either CDM subhalos with a mass Me ~ IO'^-^Mq or stars with Me ~ IMq). 
If CDM subhalos are responsible for anomalous flux ratios, then they can magnify both 
the CR and BLR because their Einstein radii (projected onto the distance of the quasar). 
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te ~ O.O1(Me/M0)^/^/;,~^/^ pc (where h = i/o/lOO km Mpc~^), are larger than the sizes 
of the CR and BLR. The NLR may be magnified as well, depending on the upper mass of 
the CDM subhalos. On the other hand, stars would be unable to magnify emission lines 
from either the BLR or NLR, while still magnifying the CR. As a consequence, depending 
on the nature of substructure, the fiuxes of emission lines in each of the lensed images that 
show anomalous (continuum) fiux ratios can be either magnified or unchanged. 

The observations of multiply lensed objects with an integral field spectrograph (IFS) 
is useful in this study, since it p rovides simultaneous information amo ng quasar images 
both spatially and spectrally (e.g.. iMetcalf et al.ll2004l : IWayth et al.ll2005l ). The IFS is more 



superior to conventional slit spectroscopic techniques in the following aspects: (1) Slit spec- 
troscopic observations would require separate observations for each pair of components. This 
not only requires more observing time but also introduces significant errors and uncertain- 
ties for relative fiux measurements since exposures for each pair are carried out in different 
observational conditions. (2) IFS allows us to optimize the relative flux measurements after 
the observations, in terms of determining the size and center positions of the apertures used 
for photometry. 

In this paper, we report on our spectroscopic mapping of the quadruply imaged quasar, 
IRXS J1131 — 1231, which shows anoma lous flux ratios, using the IFS mode of Kyoto tridi- 
mensional spectrograph II (Kyoto 3DII: ISugai et al.ll2000l . |2002| . 12004 ) mounted on the Sub- 
aru 8.2-m telescope. This lens system is nearly unique among known gravitationally lensed 
systems in that it brings together rare properties, includin g the quadruplicit y, a bright opti- 
cal Einstein ring, a small redshift, and high ampliflcation (jSluse et al.ll2003l . hereafter S03). 
S03l measured redshifts of 0.658 and 0.295 respectively for the quasar and the lensing galaxy. 
The lens shows three roughly co-linear highly magnifled images. A, B, and, C. This conflgu- 
ration emerges if the quasar is close to and inside a cusp point of the astroid caustic for an 
elliptical lens (see also Figure [7] shown in §4). In such a lens system associated with a cusp 
singu larity, there exists a uii iversal relation between the image fluxes, (Fb + Fc)/Fa = 1 
(e.g., iMao &: Schneider! Il998l ). whereas the observed flux ratios violate these rules signifl- 



cantly [{Fb + Fc)/Fa ^ 2.1 in the V band and 2.2 in the R band flS03r )]. The flux ratios 
appear to be changing with time, as a result of microlensing effects on some of the images. 
Time delay effects between the images do not play a role in the flux ratio anomaly because 
simple smooth- lens models p redict only a day or a fraction of a d ay for the time delays 
between images A, B, and C (jSluse et al.ll2006l : IClaeskens et al.ll2006[ ). Several different lens 
models with smooth mass distribution have been constructed to attempt to reproduce the 
image positions, flux ratios, and time delays of this lens system. But it appears that one 
requires a high-order lens structure in addition t o a simple ellipsoidal lens, in the forna of an 



octupole component or an unidentifled satellite (IClaeskens et al.ll2006l : [Morgan et al.l 12006 
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hereafter M06). 

Armed with the IFS mode of Kyoto 3DII, we measured the emission-hne fluxes of both 
the BLR H/3 and the NLR [OIII]AA4959,5007 for images A, B, and C simultaneously. The 
measurements of the line fluxes are more reliable than those of the continuum because 
there are no contributions from the lensing galaxy, the redshift of which differs from the 
quasar's. The H/3 and [OIII] lines are very close in wavelength, so that the effect of differential 
reddening between them is almost totally negligible. The unique fine sampling of ~ 0".l, 
together with the simple shape of the point spread function and Subaru's excellent image 
quality, provides the most suitable method for measuring the relative line fluxes. 

This paper is organized as follows. §2 describes the observations and method for data 
reduction. Spectra of images A, B, and C and flux ratios between each image pair are shown 
in S3. 54 is devoted to discussion and conclusions. 



OBSERVATIONS AND REDUCTION 



We observed IRXS Jl 131-1231 on 2005 February 8 with the IFS mode of the visitor 
instrument Kyoto 3DII, mounted on the Suba ru. We used th e atmospheric dispersion correc- 
tor installed at the Subaru Cassegrain focus (jlye et al.l 120041 ). This enabled us to obtain the 
atmospheric-di spersion-free data cube without having to apply a correction during the data 
reduction (e.g.. lSugai et al.ll2005l ). The IFS mode uses an array of 37 x 37 lenslets, enabhng 
us to obtain spectra of ~ 10^ spatial elements. The spectral range from 7300A to 9150A was 
observed in each o f two one-hour exposures. With the spatial sampling of 0".096 lenslet"^ 
( ISugai et al.l 120041 ). the field of view of ~ 3" covered the three brighter quasar images. The 
target was located at the same position of the lenslet array during the two exposures. Since 
the positional difference derived from the reconstructed data cubes was small, within a few 
tenths of a lenslet (i.e., a few xO".01), we combined them without any offset between them. 
The spatial resolution was stable during the observations, and was determined to be 0".5-0".6 
from a strong "point" -like H/5 emission-line distribution in the quasar image A (see §[3]). 

We used halogen lamp spectral frames for the flat fielding. These frames were obtained 
immediately before and after each of the target frames, with the same optical setting as the 
target frames. Since our calibration optics exactly simulates the telescope optics, i.e, since it 
provides the same micropupil images as the telescope does, it properly corrects differences in 
the lenslet-to-lenslet response (and efficiency). The uniformity of a flat-fielded reconstructed 
image is 1.5% (1 a). The wavelength calibration was carried out using He and Ne lamps, 
whose spectra were also taken immediately before and after each of the target frames. We 
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measured that the cahbration lamp hne positions changed only by 0.1 pixel (~ 0.3A) along 
the spectral direction during the whole sequence of observations. This corresponds to a 
systematic un certainty in the wa velength calibration. Moreover, the technique of micropupil 
spectroscopy (ISugai et ai.1 booeh provides us with an accuracy in the relative wavelength 
calibration of ~ 0.2A, which has been determined within each of the calibration lamp frames. 
The full width at half maximum of the instrumental line profile was measured as 7 A, which 
corresponds to 260 km s~^. In order to improve signal-to-noise ratios in target spectra, we 
smoothed them in wavelength by convolution with a gaussian profile so that the resultant 
velocity resolution was 390 km s^^. The sky subtraction was carried out using simultaneous 
sky spectra. This is particularly important for longer wavelengths, where the sky emission 
lines are strong. For an accurate sky subtraction the Kyoto 3DII has two separate fields of 
view: one for the target and the other smaller one for sky, which is located 33" away from 
the target (ISugai et al.ll2006l ). The averaged sky spectrum was subtracted from each target 
lenslet spectrum (see §[3]). 

The IFS data of the spectroscopic standard star HD93521 were used for the instrument 
response curve correction as well as for absolute flux density calibration. The extinction curve 
at Mauna Kea (CFHT Observatory Manual, http://www.cfht.hawaii.edu/Instruments/ Ob- 
servatoryManual/) was used to correct for the effects of airmass in both the target object 
and standard star spectra. The airmasses were 1.38 and 1.22 for the target frames and 1.23 
for the standard star frames. The atmospheric conditions were checked using a guide star 
count log and were determined to have been non-photometric. We roughly estimate the 
uncertainty of absolute flux density calibration to be about ±20%. 



3. SPECTRA OF QUASAR IMAGES A, B, AND C 

Figure [1] shows line images and a line-free continuum image that have been obtained 
from the IFS data. We find an extended component of the [OIII] emission connecting quasar 
images A and C. Such an extended component is not seen between quasar images A and 
B. The spatial distribution of the line-free continuum emission is similar to that of the H/5 
emission rather than the [OIII] emission. Despite the severe contamination from the host 
galaxy (§ 13. ip . the continuum is much fainter in image C than in image B. Figure [2] (upper) 
shows the spectra of quasar images A, B, and C that have each been extracted with an 
8-lenslet, i.e. a circular aperture with a diameter 0".77. An aperture of this size includes 
~ 60 % of the total flux of a point source. The aperture centers have been determined 
from the H/3 image shown in Figure [TJ The relative positions of these centers are consistent. 



to within several milliarcseconds, with peak positions obtained by lM06l using the Hubble 
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Space Telescope. As an example, the lower panel in Figure [2] shows the spectrum of image A 
both with and without sky subtraction in order to demonstrate the effectiveness of the sky 
subtraction as described in § [21 



3.1. Comparison Between A and B 

If all the emission from the quasar is magnified in the same manner within each quasar 
image, we should be able to completely subtract the spectrum of one quasar image with the 
spectrum of another quasar image by scaling the latter with an appropriate factor. If this is 
not the case, e.g. when the NLR line is magnified only by a macrolens while the BLR line is 
further magnified by a milli/micro-lens, it should be impossible to subtract these two lines 
completely with a single multiplicative factor. Using such spectral fittings, we derive the 
flux ratios in the BLR and in the NLR below. We do not attempt to derive the flux ratios 
in the continuum, which is severely contaminated by the quasar host galaxy. 

Figure [3] (upper) shows a comparison between quasar images A & B. The spectrum of 
image A has been fit, using the least squares method, in the [OIII] and the line-free continuum 
emission wavelength regions according to the following equation: /a (A) = &o x /b(A) — 6i x 
(A/8000A)^2^ where and /b(A) are the fiux densities of images A and B, respectively. 

The fitting parameter determines the scaling factor, while hi and 62 are used to remove 
the smooth continuum. The spectrum has been fit with 60 = 1-63. Similarly, we have fit 
the spectrum of image A in the broad H/5 and the continuum wavelength regions (Figure [3] 
(lower)). The broad H/5 in image A seems to have only slightly weaker emission in the bluer 
wing, compared with that in the scaled image B. The best fit 60 is 1.74. We therefore conclude 
that the fiux ratios A/B in the [OIII] and broad H/3 emission lines are 1.631qq2 ^^"^ -'-•'''^^llo 12 
respectively, where the uncertainties are derived as discussed in § 13.31 In order to compare 
these values with the ratio predicted by macrolens models, we have co nstructed a smooth- 



lens niodel with a singular isothermal ellipsoid plus an external shear (IKassiola fc Kovner 
1993 : Kormann et al. 1994), based on the quasar image positions obtained by Mod . A slight 



offset of the halo center with respect to the lensing galaxy center has been allow ed. This 



mo del predicts an A/B ra tio of 1.66. The smooth-lens models constructed by IS03l and 



by iBlackburne et al.l (120061 . hereafter BPR06), with a singular isothermal sphere plus an 
external shear, also predicts a similar ratio of 1.75 - 1.70. We conclude therefore that the 
measured ratios are close to those predicted by smooth-lens models. The required offset for 
the spectral fitting, i.e. non-zero 61, is likely caused primarily by contamination from the 
quasar host galaxy continuum. The contamination from the lensing galaxy is small since it 



is located outside of our field of view: ~ 2".l (e.g., lM06l ) down from image A in Figure [H 
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Figure [T] (right) suggests that the contribution of the lensing galaxy is less than 6 x 10^^^ 
erg cm~^ s~^ and actually much smaller since we do not detect any signature of it even 
at positions that are closer to its center. 



3.2. Comparison Between B and C 

We have compared the spectra of images B and C using the same method as in § 13.11 
When we scale the spectra according to the [OIII] line emission (Figure H] (upper)), we find 
differences in the H/3. The H/5 profile also differs in these two spectra. Figure H] (lower) 
shows the same comparison but with the C spectrum scaled according to the broad H/? 
emission. We find differences in the [OIII] lines and in the narrow line (NL) B.f3. The 
redshift of the NL R/3 matches that of the [OIII] line (z heliocentric = 0.6542 ± 0.0001 (la), 
where the uncertainty includes the accuracy of [OIII] line peak determination as well as the 
wavelength calibration uncertainties). The "NL-subtracted" spectrum in Figure H] (upper) 
suggests that the broad line (BL) H/5 is more symmetric than the original BL+NL H/3, and 
is redshifted with respect to the NL R/S {zheUocentric = 0.6581 ±0.0001 for BL H/3). We have 
measured an asymmetry parameter. A, of the H/3 emission for both the original and residual 
spectrum of image B: A = {HWHMRed - HWHMBiue)/ FWHM, where HWHMned and 
HWHMbiuc represent the half- w idths at half-max imum flux density of the red and blue 



wings of the profile respectively ( Wills et al.lll993l ). Our measurements have shown that 



A=0.32 and A=0.08 respectively, for the original BL-I-NL H/? and the residual BL H/?. The 
difference in the H/3 profile between images B and C is caused by different fractions of NL 
contribution. The above comparisons lead us to conclude that the flux ratio between images 
B and C largely depends on where the lines originate: the BLR or the NLR. The flux ratio 
C/B in the [OIII] line is 1.19+°;°^, whereas it is 0.46+°;°^ in the BL H/? line. As before, 
uncertainties are derived as in § 13.31 Compared to predictions of the smooth-lens models, 
this is larger only by a factor of ~ 1.2 in [OIII] but is smaller by a factor of ~ 2.1 in the BL 
H/3. Table [T] summarizes the flux ratios obtained between the quasar images A, B, and C, 
together with the ratios expected from the smooth-lens models. 

In order to derive the H/5 and [OIII] flux ratios between quasar images, we have used 
the simple scaling and fitting methods above, without using any template for the Fell line 
emission. This is because we have foun d that the observed Fell line ratios differ slightly from 



the well-established I Zw 1 template ( IVeron-Cetty et al.ll2004l ) and that the Fell is better 



subtracted by the scaling between quasar images themselves. 

Figure [5] demonstrates in another way the suitability of our scaling method. The H/5 
residual image has been created by subtracting three point sources with the H/5 flux ratios 
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derived above, from the H/3 image shown in Figure [H The residuals have an rms of only 3% 
of the peak value of quasar image A. 

3.3. Uncertainty Estimates 

The uncertainties of the measured flux ratios in Tabled] have been estimated as follows. 
First we have tried another spectrum fitting procedure that differs from the procedure used 
in Figures [3] and m In this new fitting procedure, we fit the line emission after removing the 
continuum emission from each image spectrum, rather than fitting the line and continuum 
emission simultaneously. This altered the flux ratios slightly, by 0.00-0.05. Secondly we have 
estimated the mutual flux contamination between images due to point-like components, 
by using the model of three point sources shown in Figure [51 In the 0".77 aperture the 
contamination in broad H/3 for images A, B, and C are 1 %, 3 %, and 5 % respectively. The 
effects on the A/B and C/B ratios are —2 % and +2 %. In the same aperture the mutual 
contaminations between the "point-like" [OIII] components (see § S]) are estimated as 2 %, 3 
%, and 2 % respectively. The effects on the A/B and C/B ratios are —1%. Thirdly we have 
derived the flux ratios with the simultaneous (line + continuum) fitting procedure, but using 
successively smaller apertures, down to the size of a 4-lenslet (0".38 diameter). While the 
C/B ratio in [OIII] decreases by up to 0.07 as the aperture size is reduced, the other ratios 
vary by no more than ±0.02. These differences should in practice include various effects: 
mutual contamination, contributions from the spatially extended components, readout and 
photon noise, and fiat fielding uncertainties. Lastly we have considered the uncertainty of the 
A/B ratio in the broad H/3 emission due to the slight difference in its wing profile between 
images A and B (§ 13. ip . When, during the fitting, we reduce the weighting of each blue 
wing data point by half, the A/B ratio is 1.76. When we reduce the weighting of each red 
wing data point by half, on the other hand, the A/B ratio is 1.69. We have taken all the 
above into consideration when determining the uncertainties of the measured flux ratios in 
Table [H The effects of uncertainties of the aperture center positions are negligible because 
the positional uncertainties are smaller than one tenth of a lenslet. 

4. DISCUSSION AND CONCLUSIONS 

Anomalous flux ratios in lensed images can be induced by substructure in a host lens. 
As a reference for investigating such substructure lensing, we estimate the angular size 
of the radius of a line-emitting region, 6s, and of an Einstein radius for a lens substruc- 
ture, 6'e, based on the set of cosmological parameters fl = 0.3, A = 0.7, and h = 0.7. 
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A source image with radius Rs at the redshift 0.658 of the quasar has an angular size of 
6*3 — 1.4 X 10~^(i?s/lpc) arcsec. For comparison, a lens having a mass Me inside an Einstein 
angle yields 9e ^ 6.7 x 10~'^ {Me/O.IMqY^^ arcsec for the lens redshi ft of 0.2 95. Thus, a 
dimensionless ratio ^ = 9s /9e sho wing the degree of a lensing effect (IBPROq ) is given as 



^ ~ 200(i?s/lpc)(ME/0.1Me 



e) 



-1/2 



Wyithe et al.l (120021 ) showed, based on detailed lensing 



simulations, that 9e must be at least an order of magnitude smaller than for no magni- 
fication. This leads us to require that ^ > 10 if there is no additional magnification of an 
image by a lens substructure. Otherwise an image should be subject to flux anomaly. 

As presented in Tabled], the flux ratios for the [OIII] line emission, which originates from 
the NLR with Rs > 100 pc, are basically in agreement with those of a standard smooth lens 
consisting of an elliptical lens and an external shear. Also, the [OIII] flux ratios appear to 
be inconsistent with lM06r s lens model, which posits a massive satellite with ~ 5 x IO^^Mq 
near image A in order to reproduce their time delay measurements. Such a massive lens 
perturber would affect the [OIII] flux significantly. The basic agreement of the [OIII] flux 
ratios with the smooth-lens model predictions suggests not only the validity of the smooth- 
lens models but also the absence of millilensing effects. A star with subsolar mass is unable 
to magnify the NLR because ^ ~ 2 x 10'^ for Rs = 100 pc and Me = 0.1 Mq. A more 
massive substructure such as a CDM subhalo, even if it exists near the light path to the 
NLR, is tightly constrained to have Me < 10^ Mq. 

It is worth noting that the observed flux ratio C/B in the [OIII] line emission slightly 
differs from the model predictions at about the 20 % level, while the flux ratio A/B is 
well reproduced. This slight inconsistency with a smooth-lens model might be caused by a 
substructure with an extended instead of p oint-like mass distribution being located at one 
of the lensed images. Ilnoue fc Chibal (120051 ) showed that for a substructure with a singular 
isothermal density distribution, the flux ratio can be changed by about 20 % from a smooth- 
lens prediction even if ^ > 10, provided the substructure is just centered at the source 
image. Alternatively, and perhaps more li kely, is the effect of an asymmetric structure or 
light distribution intrinsic to the NLR (e.g.. ISchmitt et al.ll2003l ). In contrast, a smooth-lens 
model assumes a uniform and circularly-symmetric NLR as a source image. In fact, this 
effect can partially be seen as an extension of the [OIII] emission between images A and C as 
well as a slight extension of image B in the direction opposite to images A and C (Figured]). 
This indicates that we have spatially resolved the NLR itself. Figure [6] emphasizes these 
extensions in the subtraction of three "point" source models from the [OIII] line image. 
The total [OIII] flux of the extended emission between images A and C is about a half of 
that of the point-like component of image A at our resolution. The fractional contributions 
from the extended emission components are 16 % - 20 % even in the 0".77 apertures of 
the quasar images. This causes the slight aperture dependence of [OIII] flux ratios between 
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images (particularly the C/B ratio) as described in § I3.3[ We have found, based on the lens 
modeling, that asymmetric configurations of lensed images can be caused by asymmetric 
structure in the NLR in the north-south direction, i.e. if the north side of the NLR is more 
luminous or more spatially extended than the south side. An example of such models is 
shown in Figure [71 This asymmetry in the lensed images of the NLR would be attributed to 
a slight difference between the measured flux ratios inside a flnite aperture and the smooth- 
lens predictions based on a uniform and circularly-symmetric source. 

In contrast to the [OlII] line emission, the broad-line H/5 emission, which originates from 
the BLR, shows a signiflcant difference to the model for the flux ratio C/B , i.e. C/B= 0.46 
compared with the predicted ratio of 0.91 ~ 1.00, while the flux ratio A/B is well reproduced 
by a smooth lens. This anomaly for the flux ratio C/B in the B.f3 line emission is caused 
neither by effects of (i) dust extinction nor (ii) time delay on lensed images, as argued below, 
(i) There are no effects of dust extinction on the H/? line emission, because they are not 
seen in the [0I1I]AA 4959,5007 lines, whose wavelengths are close to that of H/3. (ii) The 
time delay between lensed images B and C is a smal l fraction of a day as estimated by 



our smooth-lens model as well as by IS03I and IBPR06I . This is too short for the BLR to 
change its line flux by a factor of two, because reverberation mapping (e.g., |Petersonlll993l ) 
indicates the sizes of BLRs in AGNs to be 10^~^ light days. For these reasons, the most likely 
explanation for the anomalous flux ratio C/B in the H/3 line emission is micro/milli-lensing 
of image C. Microlensing of image C would be of particular importance, because of the finite 
probability of de-amplification even though it is located at a minimum in the arrival time 
surface (Schechter & Wambsganss 2002, see also Keeton 2003 for the case of millilensing) . 
We also note that simultaneous amplification of both images A and B by similar amounts is 
another possibility, but it ap pears to be unlikely because of the lack of such time variation 



in their continuum emission (IMOGI ). 



In addition to micro/milli-lensing of image C, microlensing of ima ge A h as occurred such 



that its CR is microlensed, while its BLR is only partially affected. IM06I have shown that 
the brightness of image A in optical bands has been increasing over recent years, probably 
because of reco very from a de- amplified stage due to microlensing. Image A was once fainter 



than image B (jS03l ) . Using Table 1 of IMOGI . we have found that the i?-band fiux ratios at 



the time of our observations (HJD=2453411) read A/B~ 1.2 and C/B~ 0.5, implying that 
the CR for image A is microlensed. As shown in § 13.11 an asymmetric profile of the H/3 
line for image A, i.e., weaker blue wing than red wing, may also indicate the presence of 
microlensing for the BLR as well but only in part; such an asymmetric profile can be caused 
by the microlensing of the part of the BLR that is rotating away from us. Thus, the Einstein 
radius of a substructure (most probably a star) near image A must be small compared with 
the size of the BLR. In contrast, a substructure located near image C would have to be more 
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massive than that near image A, because it would have to affect the emission from both the 
CR and BLR of image C, since the flux ratios C/B in both regions (~ 0.5) are almost the 
same. 

The likely radius of the BLR, -Rblr, is estimated as about 0.01 - 0.05 p c, using the 



relatio n between -Rblr and the intrinsic luminosity of H/3, L(H/5), obtained by iKaspi et al. 
(120051 ). For L(H/?), we estimate L(B.(3) ~ (1.8-4.3) x 10^^ erg s~^ as derived from the observed 



H/3 flux (which is corrected for the limited aperture size and includes the uncertainty of the 
absolute flux calibration) after taking into account (i) the dust extinction in the lensing 
galaxy and (ii) the lensing magnification predicted by lens models, as detailed below, (i) 
Our estimate of dust extinction in the le nsin g galaxy, which is an e lliptical as deduc e d from 
the presence of typical absorption lines (IS03l ). is based on work by lElfasdottir et al.l (120061 ). 
They derived, using their sample of mostly early-type lenses, that the mean differential 
extinction for the most extinguished image pair for each lens is Ay ~ 0.56 mag {Ry ~ 2.8). 
We thus adopt Ay ~ - 0.56 mag as a possible range of the absolute extinction, which 
corresponds to Ao.62/im ~ - 0.48 mag in the rest frame of the lensing galaxy, (ii) For 
the lensing magnifi cation, we adopt the predictions of smooth-lens models, yielding 12.4 
(this work) to 14.7 (jSluse et al.ll2003l ) for the magnification factor of image B. Apparently, 
L(H/5) estimated in this manner is subject to some systematics and the derived -Rblr also 
i ncludes uncertain ties associated with an intrinsic scatter in the -Rblr vs. L{}if3) relation 
( iKaspi et al.ll2005l ). Nonetheless, the range of the estimated -Rblr of IRXS J1131-1231, 
0.01 - 0.05 pc, seems to be comparable to that of a luminous Seyfert 1 rather than the 
typical size in luminous quasars (~ 0.1 pc t aken from Fig. 2 of Kaspi et al. 2005). This 
is actually consistent with the work by |S03| who, based on the estimation of its unlensed 
absolute magnitude, suggested that the source object in the IRXS Jl 131- 1231 system is a 
Seyfert 1. 

Based on our estimate of -Rblr = 0.01 - 0.05 pc, the condition ^ < 10 for the substructure 
lensing of image C's H/3 line emission requires Me ^0.1 Mq for the mass of a substructure 
near image C. On the other hand, for a substructure near image A, Me must be as small 
as 0.01 - 0.1 Mq. If the mass were larger, the H/3 flux of image A would significantly differ 
from the smooth-lens predictions, which is not the case (see Tabled]). 

In summary, we have observed the quadruply lensed quasar IRXS Jl 131-1231 with the 
IFS mode of the Kyoto 3DII. The simultaneous observation of the three brighter quasar 
images at high spatial resolution has enabled us to measure accurate flux ratios between 
them in the H/3 and [OIII] lines, and to clarify the cause of anomalous flux ratios also 
seen in the continuum emission. We have found that the flux ratios in the NLR [OIII] 
line can be explained by smooth-lens models, with only a slight deviation from the model 
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predictions. The deviation is likely caused by asymmetric structure in the NLR. The absence 
of microlcnsing and/or millilensing effects on the [OIII] line emission sets a tight limit of 
Me < IO^Mq along the light path to the NLR. The BLR H/3 shows a C/B line flux ratio 
that is smaller than the smooth-lens model predictions. This is most likely caused by the 
micro/milli-lensing of image C. The slight difference of the broad H/3 line profile in image 
A compared with those in the other images suggests the presence of a small microlcnsing 
effect on image A. Our results demonstrate that the IFS observations of gravitational lens 
systems are useful for understanding the mass distribution of lensing galaxies and the internal 
structure of lensed sources. 
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Table 1: Relative flux ratios among quasar images A, B, and C. 



Line^ 


A 


B 


c 


[OIII] 


-'-■'J '-'-0.02 


1.00 


1 iq+0.03 

-'-•-'-^-0.12 


Ep (broad) 


1 74+0.07 

' ^-0.12 


1.00 




SIEx+^ 


1.66 


1.00 


0.91 


SISx'^ 


1.75 


1.00 


1.00 


SISx'i 


1.70 


1.00 


0.96 



"Observed on 2005 February 8. Determined from the scaling and fitting procedures between quasar image 
spectra by using an 8-lenslet (0".77) diameter pure circular aperture each. See text for the uncertainties of 
flux ratios. 

''Our smooth-lcns model wit h a si ngular isothermal ellipsoid plus an external shear, based on the quasar 
image positions obtained by IMOGI . A slight offset of 9 milliarcseconds has been applied for the halo center 
with respect to the lensing galaxy center. 

■^Singular isothermal sphere plus an external shear model bv lSOSl . 
''Singular isothermal sphere plus an external shear model bv lBPROG . 
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Fig. 1. — Images extracted from our IFS data, (left) Continuum-subtracted H/3(+FeIIA4924) 
line image. The spectra have been integrated for the observed wavelength region from 7947A 
to 8I62A. The continuum level has been determined in wavelength regions from 7878A to 
7947A and from 8399A to 8469A. The length of a bar to the upper right is 1". The apertures 
used for obtaining spectra for quasar images, such as spectra in Figure [2l are shown with cir- 
cles. North is at 102° clockwise from the top. (middle) Continuum(+FeIIA5018)-subtracted 
[OIII]A5007 line image. The spectra have been integrated for the observed wavelength region 
from 8247A to 8300A. The "continuum" level has been determined in wavelength regions 
from 822lA to 8247A and from 8300A to 8343A. (right) Line-free continuum image. The 
spectra have been integrated for the observed wavelength regions from 7878A to 7947A and 
from 8399A to 8469A. 
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Fig. 2. — (upper) Spectra of quasar images A(solid line), B(dashed line), and C(dot-daslied 
line) that have been extracted with an 8-lenslet (0".77) pure circular aperture each. The 
apertures are shown in Figure [H The spectra, as well as those in Figures [3] and HI have not 
been corrected in flux density due to the finite aperture. The shorter wavelength part of 
image C is not shown, because image C was located close to an edge of the lenslet array 
and its shorter wavelength part was outside of the CCD frame when the micropupil images 
corresponding to image C were dispersed by a grism. We have carefully removed such a 
wavelength region so that the image C spectrum shown does not suffer from such flux losses 
due to the grism dispersion. Similarly, the longer wavelength part of image B is not shown. 
(lower) Spectrum of image A is shown for both cases with and without sky subtraction. 
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Fig. 3. — (upper) The spectrum of image A (thin sohd hne) has been fit with the following 
equation: = x — 6i x (A/SOOOA)*^^ where /a(A) and /b(A) arc the fiux 

densities of images A and B, respectively, and 6o, 6i, and 62 are the fitting parameters. 
The wavelength regions used for the fitting, 7880-7920A, 8160-8340A, and 8415-8455A, are 
shown with shadows. These correspond to the [OIII] and the line-free continuum emission 
wavelength regions. The dashed line is the best fit (60 = 1-63, 61 = 1.69 x 10"^^, and 
62 = 1-9), while the thick solid line shows the residuals, /^(A)— (best fit). The value of 62 is 
rather meaningless since the wavelength ranges used for the fitting are too narrow to give 
its accurate value. This causes apparent residuals in the Fell regions, which are the outsidcs 
of the fitting wavelength regions, (lower) The same as the above, except the wavelength 
regions used for the fitting and the resultant fitting parameters. The wavelength regions 
are 7880-801lA, 8071-8160A, and 8415-8455A, which correspond to the broad H/? and the 
line-free continuum emission wavelength regions. The best fit parameters are 60 = 1.74, 
61 = 2.26 X 10-^^ and 62 = -1-6. 
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Fig. 4. — (upper) The spectrum of image B (dashed hne) has been fit with the following 
equation: = cq x /c(A) — ci x (A/SOOOA)''^, where /b(A) and /c(A) are the fiux 

densities of images B and C, respectively, and Cq, Ci, and C2 are the fitting parameters. 
The wavelength regions used for the fitting, 7880-7920l, 8160-8340A, and 8415-8455A, are 
shown with shadows. These correspond to the [0111] and the line-free continuum emission 
wavelength regions. The dot-dashed line is the best fit (cq = 0.84, ci = —1.57 x 10"^'', and 
C2 = 4.2), while the thick sohd hne shows the residuals, /b(A) — (best fit). The value of C2 
is rather meaningless, (lower) The same as the above, except the wavelength regions used 
for the fitting and the resultant fitting parameters. The wavelength regions are 7880-801lA, 
8071-8160A, and 8415-8455A, which correspond to the broad H/3 and the hne-free continuum 
emission wavelength regions. The best fit parameters arc cq = 2.15, ci = 1.79 x 10""''^^, and 
C2 = —4.7. The thick solid hne shows the residuals, (best fit) — /b(A). 
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Fig. 5. — (iGft) Continuum-subtracted H/3(-|-FeIIA4924) line image. The same as that in 
Figure [H (middle) Model image with three point sources that have flux ratios shown in 
Table [1] and have the common Moffat function profile, (right) The residual image ((left)- 
(middle)). These images are shown in the same grayscale. 




Fig. 6. — Three-point-source-model subtracted [OIII]A5007 line image, where the three 
point sources have flux ratios shown in Tableland have the common Moffat function profile 
as determined from the broad H/3 point source. The apertures shown with circles are the 
same as those in Figure [H 
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Fig. 7. — The lens configuration of IRXS J1131-1231. Open circles show the observed 
positions of the images, whereas crosses show the best-fit positions of the images in our lens 
model. Solid astroid and dotted ellipse denote the caustics and critical curve, respectively. 
The center of the source position is indicated with a solid circle. In the inset, the region 
of the source position is enlarged, where a surrounding upper semicircle having a radius of 
90 pc corresponds to an asymmetric structure of the NLR in the north-south direction. This 
extended, asymmetric source image produces asymmetric arc-like features (solid curves) near 
images A, B, and C. 



